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Abstract 

Collective behavior is studied in globally coupled maps. Several coherent motions 
exist, even in fully desynchronized state. To characterize the collective behavior, we 
introduce scaling transformation of parameter, and detect the tongue-like structure 
of collective motions in parameter space. Such collective motion is supported by the 
separation of time scale, given by the self-consistent relationship between the col- 
lective motion and chaotic dynamics of each element. It is shown that the change of 
collective motion is related with the window structure of a single one-dimensional 
map. Formation and collapse of regular collective motion are understood as the 
internal bifurcation structure. Coexistence of multiple attractors with different col- 
lective behaviors is also found in fully desynchronized state. 
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1 Introduction 



Whereas the research of low dimensional chaos provided us with important 
notion of unpredictability in deterministic systems, it was soon realized that 
many natural systems are much more complicated than the low dimensional 
chaos. One of the important features in such system is high dimensionality. 
Although there remains deterministic aspects in the high dimensional chaos, 
the present nonlinear dynamics tools are not sufficient to distinguish it clearly 
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from noise. Hence, the study of high-dimensional chaos is important both from 
theoretical and practical points of views. 

Globally coupled dynamical systems, which consist of many dynamical ele- 
ments interacting all-to-all, are a good example to develop notions in high di- 
mensional systems. Such a class of dynamical systems is seen in physical, chem- 
ical and biological systems. In physics, coupled Josephson junction array[l] is a 
coupled nonlinear oscillator circuit with a global feedback. In nonlinear optics 
with multi-mode excitation [2] many modes are often coupled globally through 
energy currency. In bioscience and medical science, neural[3], cellular [4], and 
vital[5] organizations are considered as a network of active elements which are 
known to exhibit complex chaotic behaviors. Several examples in ecological 
and economic systems are also considered as a network of active agents. Glob- 
ally coupled dynamical systems is the simplest model among these complex 
network of active elements. 

So far, study of globally coupled dynamical systems has revealed novel con- 
cepts[6] such as clustering, chaotic itinerancy, and partial ordering. In partic- 
ular, study of collective dynamics has gathered much attention [7-21]. When 
the interactions between elements are small enough, each element oscillates 
independently without synchronization between them. Therefore the degrees 
of freedom of the system are effectively proportional to the system size. If each 
element has chaotic dynamics, the system may be thought as high dimensional 
chaotic state. Even in such a case, a macroscopic variable show some kind of 
complicated dynamics rather than noise, ranging from low-dimensional torus 
to high-dimensional chaos[18,12]. This may imply that any weak interaction 
between active elements necessarily brings some sort of correlation among 
them. 

The purpose of the present paper is to study the nature of such collective 
motion adopting a globally coupled map [6], and present a mechanism for the 
origin of such collective dynamics. With the change of the control parameters, 
collective dynamics shows some sort of bifurcation. We present how elements 
are organized to show the bifurcation structure in the collective dynamics. 

In Section 2, globally coupled logistic map is introduced and its characteristic 
phenomena are presented as a brief review. In Section 3, an overview of differ- 
ent kind of collective dynamics in the desynchronized state of globally coupled 
logistic maps is presented. In macroscopic dynamics, lower dimensional motion 
and much longer time scale than that of microscopic dynamics are observed. 

Our interest is focused on the thermodynamic limit of such collective behav- 
ior. In Section 4, the time scale and the amplitude of collective motion are 
studied in the limit of large system size. In Section 5, global phase diagram in 
the parameter space is presented. While the phase diagram shows a compli- 
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cated structure, tongue-like bifurcation structures are clarified by introducing 
a scaled nonlinearity parameter. Collective dynamics with a larger amplitude 
exists in each tongue structure that corresponds to a periodic window in the 
single logistic map. The elements are accumulated to few bands corresponding 
to the window for a small coupling. Since windows exist in any neighborhood 
in the parameter space, the clarification of the collective dynamics with such 
bands is necessary to understand the collective dynamics in general. Thus we 
focus on such tongue structures in Section 6, to reveal a mechanism of col- 
lective dynamics, where internal bifurcation of elements plays a key role. In 
Section 7, bifurcation of tongue structure is studied in connection with the 
internal bifurcation. Even within the same tongue structure, we can observe 
different types of collective motion. The growth of tongue structure with the 
coupling strength is also discussed. In Section 8, hysteresis and multiple at- 
tractor phenomena of the collective motion are reported. This paper concludes 
in Section 9 with summary and discussion. 



2 A Simple Network Model of Chaotic Elements on Globally Cou- 
pled Map 

In the present paper the following Globally Coupled Map(GCM) is studied. 



= (1 - e)f{xn{t)) + i- ^ f{xn{j)), (^ = 1, 2, 3, • • • iV), (1) 



where Xn{i) is the variable of the ith element at discrete time step n, and f{x) 
is the internal dynamics of each element. For the internal dynamics we choose 
the logistic map 



where a is the nonlinearity parameter. The logistic map has been studied 

in detail as a typical of dissipative chaos. The parameter e gives a coupling 
strength between elements. The total number of elements denoted as A^. The 
nonlinearity parameter a, the coupling strength e, and the system size N are 
the control parameters of the GCM. 

The GCM can be considered to be a mean-field extension of coupled map 
lattice(CML)[22], in which elements are located at discrete spatial coordinates 
and interact with neighbors. GCM can be also considered as a CML in which 
the spatial dimension goes to infinity. 
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In the GCM model, two opposite tendencies coexist: all-to-all coupling tends 
to synchronize elements, while chaotic instability in each element tends to 
desynchronize them. Depending on the balance between the two tendencies, 
a rich variety of phenomena has been found [6]. When the coupling strength 
is strong enough, all elements are synchronized each other and the dynamics 
is nothing more than the single logistic map as is called coherent phase. As 
the coupling strength is smaller or the nonlinearity larger, elements split into 
some groups, in each of which they are synchronized each other. This regime 
is called ordered phase, while the phenomena are called clustering. The 
clustering is common characteristics in globally coupled systems, including 
globally coupled oscillator systems[23,24]. 

In the region (called partially ordered phase) where the two opposite ten- 
dencies are somewhat balanced, some part of the elements makes a few clus- 
ters, while the rest elements do not form clusters and their oscillations are 
desynchronized. In the phase space, there are a lot of "attractor ruins" with 
lower dimensionality, at which the trajectory is attracted and stays over some 
duration, but then the trajectory goes out from them into much higher di- 
mensional phase space, till they are again attracted to another attractor ruin. 
In this phenomenon, called chaotic itinerancy, effective degrees of freedom 
changes with time[6,25,26]. 

If the coupling strength e is small enough, desynchronizing tendency is so 
dominant that elements are mutually desynchronized [7, 8] (called desynchro- 
nized phase). In this case, in general, all Lyapunov exponents are positive[^ 
and the degrees of freedom in this system are proportional to the system size 
A^. In the desynchronized state, the mean-field is not stationary, where a kind 
of collective dynamics is observed. In this paper we focus on the macroscopic 
state on this desynchronized state. To see the macroscopic state, the dynamics 
and statistics of the mean-field 

1 ^ 

hn = ^j:fMj)), (3) 



are studied as an order parameter. In the next section, we will show some 
phenomena of macroscopic dynamics in the desynchronized state. 
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Fig. 1. Time series and return map. Time series are plotted at every 2 steps after 
transients are discarded. The parameters are (a)a = 1.5449205, e = 0.0005, = 10^, 
and (b)a = 1.5449205,6 = 0.0005, N = W. Corresponding return maps {hn, /in+io) 
arc plotted over 50000 steps after transients arc discarded. 

3 Phenomenology of Collective Motion in Desynchronized State 



Since in the desynchronized state there is no mutual synchronization in ele- 
ments, one might imagine that the mean-field would be effectively the same 
as noise and therefore the mean-field goes to a constant with the increase 
of N. One might consider that such high dimensional dynamics can not be 
distinguishable from noise. 

Indeed this is not the case. One of the authors has found that the mean 
field dynamics is different from noise, and studied its nature as "hidden co- 
herence" [7,8]. A simple solution to such collective dynamics is the possibility 
that it is represented by low-dimensional dynamics in the thermodynamic limit 
(A^ — > oo), even though each element is chaotic and mutually desynchronized. 



^ Number of negative Lyapunov exponents may be related to the number of 
synchronization [29] . 
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Fig. 2. Power spectrum of time series of a single element (upper) and the mean 

field /i„ (lower), a = 1.5449205, e = 0.0005, N = lO'^, While the spectrum for a single 
element has a peak at the frequency 0.5, the dynamics is more irregular than the 
mean field. The slow dynamics of the mean field is shown at the frequency 0.025269. 

Indeed such examples have recently been found in short-ranged coupled map 
lattice and cellular automata[10], globally coupled oscillators [2 7, 28], globally 
coupled tent map[12, 16,20,21], and globally coupled logistic map with hetero- 
geneous elements[18]. 

In the present case, the collective dynamics is not given by such low dimen- 
sional dynamics [8], although it has some structure distinguishable from noise. 
Let us give a few sets of examples of the mean-field dynamics. 

Fig. 1(a) shows the time series of the mean- field as a function of time step 
n at every 2 steps, and the corresponding return map of the mean field (for 
a = 1.5449205, e = 0.0005, = 10^). The coupling strength is too small 
to synchronize any two elements. The trajectory of the mean-field has some 
fluctuation due to the finite system size. With the increase of the system size N 
(parameters are a = 1.5449205, e = 0.0005, = 10^), however, the trajectory 
shows some coherent motion as is shown in Fig. 1(b). The trajectory is rather 
close to quasiperiodic motion, although the points are scattered around the 
"torus" motion. In Fig. 2. power spectra for the time series of an element and 
the mean-field are overwritten. The mean-field dynamics has a much longer 
time scale than that of an element. 

Note that the width around the closed curve remains finite with the further 

increase of A^. The collective dynamics is not on a two-dimensional torus, and 
indeed is not represented by low- dimensional dynamics as will be demonstrated 
in the next section. On the other hand, since the mean-field dynamics does not 
approach a point with the increase of N, it is also different from noise. Hence 
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Fig. 3. Time series and return map. Time series are plotted at every 7 steps. The 
parameters are a = 1.69620, e = 0.008, N = lO'^(a), a = 1.69755, e = 0.008, 
N = 10^ (b), a = 1.69844, e = 0.008, = lO'^(c). See also Fig.4 

the collective motion has some structure, although it is high-dimensional. 

Another set of examples is given in Fig. 3, which are the time series plotted 
at every 7 steps and the first return maps. In Fig. 3(a) (parameters are a — 
1.69620, e = 0.008, = 10^), quasi-periodic-hke motion is not detected in the 
mean field dynamics, hut some structure exists in the return map, while in 
the time series, characteristic time scale seems to exist. 

By slight increase of a (i.e., with the parameters a — 1.69755, e = 0.008, N — 
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Fig. 4. Power spectrum of time series of a single element a;„(i) (upper) and the mean 
field /i„ (lower), corresponding to Fig.3. 
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10^), the dynamics of the mean field is changed as in Fig.3(b). In this case, 
the return map does not show a clear structure, and the variation of the mean 
field remains at the same magnitude with the further increase of A^. 

With much slighter increase of a (parameters are a = 1.69844, e = 0.008, 
= 10''), the mean field comes to oscillate more regularly, whereas the motion 
is scattered around torus motion(Fig.3). 

Note that the choice of every 2 or 7 step in the above plots is not arbitrary 
but there is a reason for it, as will be clear in the following sections. Our 
goal in this paper is to give a consistent explanation for the above collective 
motion, and answer the remaining questions in the collective dynamics: When 
the system size goes to infinity, i.e., in the thermodynamic limit, how is the 
macroscopic dynamics characterized? How does the remnant order in high- 
dimensional collective dynamics emerge out of the complete desynchronized 
elements? How is a longer-time scale in the collective dynamics formed? How 
does the collective dynamics depend on the parameters a and e, and what 
kind of bifurcation structure is expected, and how is it explained in terms of 
dynamical systems theory? 



4 Thermodynamic Limit of Collective Motions 

4-1 Amplitude of Collective Motion 

In the previous section, we have mentioned that the mean field dynamics plot- 
ted in the return map shows some structure, ranging from lower-dimensional 
structure, such as torus, to higher- dimensional stochastic structure. First of 
all, to characterize the mean field dynamics, we measure the mean square 
deviation(MSD) of the mean field distribution, 

{{6hr) = {h') - {h)\ (4) 

as a measure of the amplitude of the mean field dynamics. The bracket (■) 
denotes the temporal average. Since the motion of the mean field is not on a 
torus, it is not always possible to define the amplitude of oscillation by the 
radius of the torus pattern. Even though such collective oscillation is hardly 
detected, the above MSD works as a measure for the amplitude, and also is 
useful to measure the variation around the fixed point[^. 

Since there is no synchronization each other, the MSD also provides a simple tool 
to see whether the population obeys the law of large numbers. 
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Fig. 5. Mean square deviation (MSD) of the mean field distribution are plotted as 
a function of the system size N. 

In Fig. 5 the MSD of the mean-field is plotted with the system size N. In some 
case, the MSD decreases up to a certain constant with the increase of N, but 
remains constant with the further increase. In other cases, the increase of MSD 
is seen at some range of size, but then approaches a certain constant. These 
show the distinction of the mean field dynamics from pure noise and suggest 
some coherence between elements. 

4-2 Degrees of Freedom of Collective Motion 

Next, we study the degrees of freedom in the collective dynamics in the thermo- 
dynamic limit. In the previous section, we have mentioned that the collective 
motion, detected in the return map, has some low- dimensional-like structure 
but the width of scattered points around the 'torus' remains finite in the 
thermodynamic limit. Since the possibility of higher dimensional torus is not 
excluded only by the figure, we measure the correlation dimension[30] of the 
mean-field time series. 

In Fig. 6, the change of slope in the correlation integral ^^^f^^J^^ is plotted as 
a function of the scale size with increasing the system size N. For a smaller 
system size, the correlation dimension is increased monotonically with the 
decrease of the scale as in the case of random variable. For a lager system 
size, curves have a plateau at a value less than the correlation dimension two, 
which seems to correspond to the collective motion. In a smaller scale, however. 
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Fig. 6. Correlation dimensions are plotted as a function of scale size r for different 
system size, which are indicated at the right of each figure. The parameters are 
a = 1.699, e = 0.008(a), and a = 1.5449205, e = 0.00050(b). 
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Fig. 7. Correlation dimension for heterogeneous GCM[18]. The parameters are 
a = 1.9, e = 0.11. Nonhnearity parameter for each element is homogeneously dis- 
tributed over a £ [1.875, 1.925]. 



correlation dimension becomes large. At this smaller scale, the motion is hard 
to be distinguishable from noise. If the scale of this regime got smaller with 
the size N, one could conclude that the collective dynamics is low-dimensional 
in the thermodynamic limit. As shown in Fig.6, this is not the case. The slope 
function converges to a certain curve with the increase of size N where the 
plateau region is no more expanded. Thus, the mean field dynamics does not 
converge to lower dimensional dynamics in the thermodynamic limit. 

To check the validity of this method, it will be relevant to mention the case 
with a heterogeneous system[18], e.g., a globally coupled map with distributed 
nonlinearity parameter a over elements, where the mean-field dynamics shows 
a clear quasi-periodic motion. The width of scattered points around the tours 
converges to in the thermodynamic limit. Corresponding plots of slopes are 
given in Fig. 7, where the plateau at the value 1 is expanded with A^, and the 
"noise" region is shrieked to the scale r ^ 0. The difference from our uniform 
case is clearly visibleQ. 



° In the heterogeneous case, the law of large numbers might be considered to be 
recovered around the torus motion. 
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Fig. 8. Rotation number of the mean field dynamics R, plotted as a function of 
system size N. 



4-3 Characteristic Time Scale of Collective Motion 



To see how the time scale of the collective motion depends on the system size 
N, we have measured the rotation number of the mean field dynamics as a 
function of the system size N. Here, the rotation number R is defined as 

= lim - V ^ (5) 



where A^„ is angle variable formed by two vector — {h),hn+i — (h)), 
{hn+i — {h), hn+2 ~ {h)) defined around the average mean filed (h) over time. 

In Fig. 8, the rotation number converges to a certain value. This implies the 
appearance of characteristic time scale in the mean field dynamics, indepen- 
dently of the system size for large enough N. 

Power spectrum of the mean field dynamics also supports the existence of 
characteristic time scale of the mean field dynamics in the thermodynamic 
limit as is shown in Fig.9. While the spectra indicate that the mean field 
dynamics is non-periodic, there are peaks, which get shaper with the increase 
of N up to certain size and converge to a certain curve. 
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Fig. 10. Mean square deviation (MSD) of the mean field dynamics is plotted as a 
function of a. e = 0.02(o), 0.15(+), 0.01(n).A^ = 2^^. 

5 Global Phase Diagram of Collective Motion in Parameter Space; 
Tongue-Like Bifurcation Structures 

As we have shown in the previous section(§3), the mean-field dynamics de- 
pends on the parameter a, and e. In this section, we study how the collective 
motion depends on the parameters. To characterize this dependence on the 
parameters it is often convenient to use the the mean square deviation(MSD) 
of the mean field distribution as we have introduced in §4. 

5.1 Phase Diagram in (a, e)-Plain 

In Fig. 10, the MSDs are plotted as functions of the parameter a for several 
coupling strengths e. Here the system size is chosen to be large enough, to 
see the behavior of MSD converged in the thermodynamic limit. Two points 
should be noted here. First, the change of MSD is not monotonic with a, but 
is rather complicated. Second, the change of MSD is comphcated with fine 
structures, which still keep some similarity against the changes of the coupling 
strength e. For example, a similar but slightly different structure is visible for 
a ^ 1.7025 for e = 0.01, a ^ 1.725 for e = 0.015, and a ^ 1.73 for e = 0.02. 
In Fig. 11 the parameter dependence of MSD is plotted on the 2-dimensional 
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Fig. 11. Mean square deviation (MSD) of the mean field dynamics is plotted in (a, e) 
plane with gray scale. The scale shows the value of MSD, where the darkest one 
corresponds to MSDfsi 10"^, and the brightest one to MSDft! 10"^. 
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Fig. 12. Mean square deviation (MSD) of the mean field dynamics hn are plotted 
as functions of the effective nonlinearity parameter A. 

(a, e) plane. First, regimes with a larger collective motion form tongue-like 
structures, each of which starts at some point or intervals of parameter a at 
e = 0, and grows with e. Second, the growth of the edge in a tongue-hke 
structure has a nonlinear dependence on the parameters a and e. Third, for 
almost all parameter values, the MSD of the mean-field remains finite in the 
thermodynamic limit. 

5.2 Effective Nonlinearity 

To see the structure in the parameter space closely, we introduce rescaling of 
the parameters. For it, we note that each element obeys the following dynam- 
ics, 

Xn+i = {I - e){l - axl^) + ehn (6) 

where hn is the mean-filed value at time step n, which can be considered 
as a time dependent input to each element. In this map, the nonlinearity is 
modified by the additional term hn- Taking into account of this point, we 
normalize the variable follows: 

Xn ^ {1 - e + ehn)Xn. (7) 
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Fig. 13. (a)Mean Square Deviation (MSD) of the mean field dynamics is plotted as 
a function of the effective nonlinearity parameter A and e. The scale shows the value 
of MSD, where the darkest one corresponds to MSDsi 10~^, and the brightest one 
to MSD~ 10~^ (see Fig.ll).(b)Logistic map bifurcation diagram with the increase 
of nonlinearity parameter. Horizontal axises are common among two figures. 
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Then the logistic map of each element is written as 

= l-a(l-e)(l-e + e/i„)X„' (8) 

where a(l — e) (1 — e + e/i„) can be regarded as the nonlinearity at time step n. 

Since the above scaling is time-dependent due to hn, we define the effective 
nonlinearity parameter A as 

A = {1 - e){l - € + €{h))a. (9) 

with the time independent rescaling of 

Xn^{l-e + e{h))Xn. (10) 

where {h) is mean-field average in time, 

(M^fcig/i^. (11) 

In Fig. 12 we have plotted MSD by adopting the effective nonlinearity param- 
eter A instead of a. In Fig. 13 (a) the parameter dependence of MSD is plotted 
on the 2-dimensional (^4, e)-plane. The scahng structure of tongues seems to 
be much clearer. While the width of each tongue seems to increase roughly 
linearly with a and e, detailed discussion will be appeared in §7.3. 

When the coupling strength e approaches 0, each tongue structure corre- 
sponds to a window of the single logistic map(Fig. 13(b)). For instance, be- 
tween A ^ 1.75 and A ^ 1.79 a tongue structure is clearly shown in Fig. 13 (a), 
corresponding to the period-3 window of the single logistic map. Although 
there are countably infinite windows in the parameter space in the logistic 
map, it is difficult to detect the windows for a longer period numerically. 
However, it is remarkable that a lot of tongue structures are visible in our 
model, corresponding to the windows with a longer period. 

In each tongue structure, further internal structures exist. For instance, the 
tongue corresponding to period-3 window of the logistic map between A — 1.75 
and A ps: 1.79, has three internal structures, roughly speaking. To understand 
each inner structure in the tongue, we will study the dynamics of each element 
and the distribution in the following sections (§6, and §7). 
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Fig. 14. The distribution dynamics is plotted as a function of time. The density of 
the population is shown with the use of a gray scale. The darker region indicates 
the high density. The parameters are a = 1.8445,e = 0.015, N = 10^. 

6 Collective Behavior Through Self-Consistent Dynamics 

In this section, we briefly describe how the collective motion is formed, espe- 
cially focusing on the tongue structure. 

6.1 Distribution Dynamics 

In the limit of N oo, the probability distribution function is defined as 
follows, 

= Jim^-^^5(x-x„(i)). (12) 

i 

Oscillation, rather than the fixed point, of the mean field dynamics implies 
that the probability distribution function does not also remain stationary but 
depends on time. 

Time series of the probability distribution function by numerical calculation is 
given in Fig. 14. The parameters for the figure (a = 1.8445, e = 0.015) belong 
to the tongue structure in the period 3 window. In this case, since the mean 
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field dynamics has the component of period 3, we plot the figTire by every 3 
steps to see the slow modulation of Pn{x). Due to the chaotic oscillation of 
each element, the distribution function spreads over x G [—0.8,1.0], but the 
distribution is not monotonous, and has some structure. In the three regions 
around x 1.0,0.0,-0.8, the population are relatively large. This number 
"three" corresponds to the period of window which the tongue structure of 
this motion corresponds to. The number of elements in each three region 
oscillates in time, and the phase of each oscillation is different from that of 
the others. 



6.2 Formation of Self- Consistent Dynamics 

It is interesting to view our collective dynamics as an interference between 
mean field dynamics and individual elements. Before we present a scenario for 
collective motion, we demonstrate the formation of self-consistent dynamics 
between the mean field dynamics and individual elements as follows. 

For simplicity, we adopt the case, in which the effective nonlinearity parameter 
A is near the period-2 band merging point(the parameter are a = 1.5449205, 
e = 0.0005, and the time series and the return map for the parameter are given 
in Fig.l). The distribution function is given in Fig. 15 at every 20 steps. In this 
case, distribution of elements can be divided into two regions around x 0.54. 
During these 40 time steps the value of distribution function at the left region 
{x < 0.54), given in Fig. 15(b), decreases with time, while the other region 
plotted in Fig. 15(c) increases with time. Although the change of distribution 
is quite small, there is a systematic oscillation in the distribution(cf.Fig.20). 

Consider a population dynamics of each of the two regions. In Fig. 16(a), the 
population in each region, Nl, and Nr are plotted as a function of time. Nl 
denotes the number of elements in the region smaller than x = 0.54 in Fig. 15, 
while Nji(— N — Nl) denotes that for larger than x = 0.54. (The definition 
for each region is described below in detail.) The population in each region 
oscillates in time. In Fig. 16(b), on the other hand, since the mean field has a 
component of period two, the time series hn-i and hn are plotted at every two 
steps. Note that the mean field also oscillates in time with the same period as 
Nji, and Nl, but the phase of the mean field oscillation is different from that 
of the population dynamics in Fig. 16 (a). 

To see how the mean filed dynamics and the distribution dynamics interfere 
each other, we construct a return map of above two quantities. Fig. 17 is a 
return map of the distribution dynamics and the mean field dynamics. This 
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Fig. 15. The distribution functions at time n = 0(*), 20(«), 40(o) are shown, 
figure implies that a self-consistent dynamics is formed as follows, 



hn = h{hn-l,Nn-l), 



(13) 



where each h and is a function of /i„ = hn-i — hn, and Nn = Ni — Nf>. If the 
mean field were an external force for each element, the population responds to 
the mean field value. On the other hand, since population organizes the mean 
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Fig. 16. (a)Time series of the number of elements in the 2 regions, (b)time series of 
the mean field.a = 1.5449205, e = 0.0005, N = 10^. 
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Fig. 17. Return map of the time series of number of elements in the 2 regions 
and time scries of the mean field, plotted per two steps for even n and odd n. 
a = 1.5449205, e = 0.0005, = 10^. 
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Fig. 18. Behavior of the third iterate of a logistic map. (a) After tangent bifur- 
cation at three points, and (b)after crisis. Sohd hnes with arrow indicate iterated 
trajectories starting from three points 

filed dynamics, the collective motion can be described as a self-consistent 
relation between the population dynamics and the mean field dynamics. 

Prom the above viewpoint, we now demonstrate how the population distribu- 
tion is modified as the mean field varies slowly. 



6.3 Internal Bifurcation in Temporal Domain 



If the mean field were an external force for each element, it could be regarded 
that each element follows the logistic map with external force. This is valid if 
the mean field varies slowly. In this case, the motion of equation for the i'th 
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Fig. 19. Behavior of the third iterate of the map Eq.(15). In contrast with the case of 
Fig. 18, only one or two regions can attract trajectories. (a)The region around x = 
attracts orbits as a region after tangent bifurcation, while the other regions, which 
are before the tangent bifurcation, cannot attract orbits. (b)Two regions around 
X = 0, and x = —0.8 attract orbits, while the region around x = 1.0 can not attract 
orbits due to the crisis. (c)The region around x = attracts orbits, while the other 
regions can not attract orbits due to the the crisis. 
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Fig. 20. Temporal change of 5pn{x) = Pn{x) — {p{x)) is plotted by the thick line 
by magnifying the scale 20 times, while 6peq-n{x) = peq,n{x) — {p{x)) is plotted 
by the thin line, where {p{x)) = limj^oo j Z)n=i Pn(a^)- We have plotted them 
instead of Pn{x) and p^eq,n{x) whose change is too small to be visible. Roughly 
speaking, the population in the region where 5p^eq,n{x) is negative is going to the 
region where 5p^eq,n{x) is positive. 5p^eq,n{x) changes qualitatively at n = 12 (c). 
a = 1.5449176, e = 0.0005, N = 10^. 
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element is given by 



Xn+lii) = Fn{Xn{i)), (14) 

with 

F„(x) = (1 - e)(l - ax^) + e{h) + e ■ Sh^, (15) 
where Shn is difference from (h), i.e., /i^ = (h) + 5hn- 

If we took 6hn out of consideration, the dynamics of each element would be 
same as the logistic map with the nonlinearity parameter A = a{l — e){l — 
e + e{h)). Since in the previous section (§5) each tongue structure has good 
correspondence with a window of the logistic map, here we especially focus on 
the window structure of the logistic map. In the case of the logistic map, the 
period p window starts at the tangent bifurcation point of the p'th iterate of 
the map, and then the period doubhng bifurcation proceeds with the increase 
of a, until the window ends up by crisis(see Fig. 18). Note that for a period p 
window tangent bifurcation or crisis occurs at p points of x at the same value 
of A. In this case, since Eq.(15) has an invariant measure, the probability 
distribution in each of p pieces of regions is equivalent. 

Take 5hn in Eq.(15) into account as an external force. The bifurcation of a 
logistic map with time dependent external force has a crucial difference from 
usual bifurcation of the logistic map. In Fig. 19, examples of the third iterates 
of the map with external force are shown. In Fig. 19 (a) a region around x ^ 
crosses y = x, while two regions around x ~ 0.95 and x ~ —0.75 do not cross 
y — X. In Fig.l9(b)(c), while three regions cross y — x, one or two of the 
regions are collapsed by crisis. 

In general, consider the case of period p window with external force. Since the 
tangent bifurcation or crisis of p points occurs at a different value of A, the 
number of divided attractors can be changed with A. Even if trajectories are 
attracted into p distinct regions, the probability distribution in each region is 
not equivalent. 

As we have seen in the previous subsection (§6. 2), slow modulation of the mean 
field leads to the dynamics of the distribution of population. With the slow 
modulation of Shn in time, behavior of each element also changes. In other 
words, with the change of Shn, bifurcation can occur in the effective map for 
each element, 

= Fn+p-l O Fn+p-2 o • • • O o F„, (16) 
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which is the p'th iterate of the map Eq.(15). Since Shn changes temporally, 
such bifurcation occurs temporally. To distinguish from the notion of bifurca- 
tion in parameter space, such bifurcation is called as "internal bifurcation" 
While the notion of "internal bifurcation" indicates the slow modulation of the 
effective map Eq.(16), since time dependence of 6h induces asymmetry in this 
effective map, as we have shown in the previous paragraph, the dynamics of 
the mean field has component of period p motion. This is why we have plotted 
the time series of the mean field and the probability distribution function at 
every p steps. 

To characterize the effective map at every p time steps, we introduce the 
invariant measure p^fjni^) determined from the Perron- Frobenius equation, 

1 

pit(^) = / pSn(?/)'^(^ - Fjr\y)))dy- (17) 

If the mean field changes slowly, we can approximately regard our GCM dy- 
namics as relaxation process of Pnix) to Pegln(^)- each, time step, elements 
follow the effective map Eq.(16) so that the distribution function /9„(x) is go- 
ing to be relaxed. The elements in the region where p^fgni^) — pn{x) < are 
going to the region where Pegln(^) ~ Pn{x) > 0. As a result pn{x) is going to 
relax into p'f^ni^)- 

On the other hand, the mean field is derived as, hn = J f{x)pn{x)dx. Relax- 
ation of pn{x) can lead the mean field to a certain critical value, at which the 
internal bifurcation occurs in the effective map Eq.(16). For instance, small 
difference of the mean field induces one point in the effective map to be tan- 
gential toy = X, or one region to be collapsed by crisis. As a result, the nature 
of the invariant measure p^fjni^) effective map changes qualitatively. 

With this internal bifurcation, the distribution Pn{x) is not actually relaxed 
to p^Jg'nix), because 1) the velocity of change Pn{x) is finite, and 2) the relax- 
ation of Pn{x) makes Pegln(^) t)e changed qualitatively. Consequently Pnix) 
oscillates in time. This is a qualitative explanation why the mean field does 
not approach a fixed pointQat the thermodynamic limit. 

^ In our previous work [18], the nonlinearity parameter a was distributed over 
elements. In that case, some sort of differentiation of dynamics over elements enabled 
the collective motion possible. To characterize the differentiation, the notion of 
"internal bifurcation" was introduced as a snapshot representations of one system. 
As we will show below, since the temporal bifurcation in an element leads to the 
collective motion, we extend the notion of "internal bifurcation" to identical case. 
^ The unstable fixed point of the mean field value is given as h = J f{x)po{x)dx, 
where po{x) is a fixed point solution of Eq.(17) with p = 1. The fixed point solution 
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Let us look back to the example in the subsection §6.2 and try to describe the 
dynamics along the above scenario. Population distribution function p„(x) at 
time step n is plotted with the solid line in Fig.20(and see the caption in it). 
Since the efTective nonlinearity parameter A is near the band merging point of 
the logistic map, it is useful to define the two regions as follows. The effective 
map given by the second iterate of map, 

Fn{Fn-i{x)) = (1 - e)/((l - e)f{x) + e/i„_i) + eK, (18) 

has three unstable fixed points and the middle of these points is denoted by . 
R and L denote the region where a; > x* and x < x* respectively. (Based on 
this definition, we have calculated the number of elements in the two regions, 
from which Fig. 17 in the subsection §6.2 is obtained.) 

In this parameter, if the unstable fixed point of the mean field solution were 
realized, these two regions collapse due to the crisis. Hence, these two region 
are unstable. (In the bellow "stable" or "unstable" means that a region, R 
or L, can attract trajectories or not.) As we have discussed above, however, 
dynamics of the mean field modulates the effective map Eq.(18), and then, 
for this parameter regime, there are three cases: 1) R region is stable and L is 
unstable, 2) R region is unstable and L is stable, or 3) both L and R regions 
are unstable. 

The effective map Eq.(18) can be characterized as invariant measure of ef- 
fective map Peq-'„(x), which depends on the mean field value. The stability in 
each region can be seen by the strength of p^^Jni^)- The thin lines in Fig. 20 
are p^^J^{x){see the caption in Fig.20). In Fig. 20(a), R region is unstable and 
L region is stable at n = 1 as is shown by piq^„(x). Since p„(x) is going to 
relax to piq|„(x), the elements in R region moves to L region as is shown in 
Fig. 20(a) and (b). Indeed the number of elements decreases in R region and 
increases in L region in Fig. 16(a) and (b). This change of p„(x) continues 
until the modulation of the mean field induces crisis of L region at n = 12 in 
Fig. 20(c). By this destabilization of L region, elements in L region move to 
R region so as to relax the distribution pn{x) to piq^„(x) in Fig. 20(d) and (e), 
until the next crisis leads to the structure of Fig.20(a), giving a flow from it! 
region to L. 

To sum up, distribution function p„(x) changes slowly so as to relax to peq,n{x) 
and the distribution Pnix) changes until the modulation of the mean field in- 
duces internal bifurcation structure to change qualitatively. In this example, 
qualitative change in internal bifurcation is due to local crisis. By the change of 
stability in two regions, p„(x) relaxes to a different region. With the repetition 



Po{x) is unstable. 
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Fig. 21. Ratio of the time steps, in which the two regions are both sta- 
ble("stable+stable"), one of the regions are stable("unstable+stable"), and the two 
regions are both unstable("unstable+unstable"). e = 0.0005, N = W . The insta- 
bihty is due to the crisis bifurcation in internal bifurcation. 

of this stability change the mean field oscillates quasi-periodically. This mecha- 
nism of the stability change in each band holds for any period-p band (window) 
regime where elements are attracted to and repelled from each band region 
successively with the internal bifurcation giving rise to crisis. 



7 Bifurcation of Collective Motion 

7.1 Bifurcation of Tongue Structures 

As wc have seen in Fig. 20 in the previous section (§6.3), one of the two regions 
in the second iterate of the effective map Eq.(18) is collapsed due to the crisis 
at some time steps and such a region changes in time. With the increase of 
a, the time interval of crisis bifurcation becomes longer. In Fig. 21, the ratio 
of the time interval, during which one of the two regions is unstable and the 
other is stable, the two regions are both unstable, and the two regions are 
both stable, are plotted with the change of the nonlinearity parameter a. For 
a < 1.5449173, crisis bifurcation never occurs both in the two regions, while for 
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Fig. 22. Mean square deviation (MSD) of tlie mean field distribution with the in- 
crease of e plotted as a function of the effective nonlinearity parameter A. The region 
corresponds to the period-2 tongue structure. Aq denotes the band merging point 
(crisis bifurcation point of 2-band) of the logistic map (^o = 1.543689012692076). 

a > 1.5449173, the time interval of crisis is getting longer. For the parameter 
beyond a = 1.5449212, the two regions is collapsed due to the crisis bifurcation 
all time steps. Hence, the period-2 tongue structure starts at the parameter, 
where one of the two regions in the second iterate of the effective map Eq.(18) 
is collapsed due to the crisis at some time steps at a = 1.5449173, and ends 
up at the parameter where both the two regions collapse due to the crisis all 
the time at a = 1.5449212. 

Consider an internal bifurcation condition of Eq.(16) (for instance, crisis bifur- 
cation or tangent bifurcation in each element.). While for e = the bifurcation 
condition holds at only one point in the nonlinearity parameter, for finite e, 
due to the oscillation of the mean field, the internal bifurcation condition is 
satisfied for some steps within some interval Asmaiii^) < ^(e) < Aiargd^) in 
the parameter space. Hence, the edge of a tongue structure, corresponding to 
a periodic window of logistic map, starts from tangent bifurcation and crisis 
bifurcation point at e = 0, and each line constitutes the parameter A and e, 
where each bifurcation condition holds at some time steps(Fig.22 for period- 
2 tongue structure). Scaling the width of tongue structure will be discussed 
lated (in §7.3). 
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1.2 Bifurcation in a Tongue Structure 

Even within the same tongue structure, we can observe different types of col- 
lective motion. With the change of the parameter a and e, the collective dy- 
namics shows a kind of bifurcation. Since the collective dynamics remains high- 
dimensional, it is not described as a standard bifurcation in low-dimensional 
dynamical systems. Here we study a mechanism of such change in the collec- 
tive dynamics. 

In Section 3, it is shown that slight increase in a induces the qualitative change 
of collective dynamics (Fig. 3). To see this quantitatively, it may be convenient 
to measure the rotation number of collective dynamics. In Fig. 23, the rota- 
tion number is plotted as a function of a. In the regime plotted in the figure 
(i.e., between a G [1.69614,1.69847] with e = 0.008), period-seven tongue 
structure is observed. Roughly speaking this period-seven tongue region is di- 
vided into 3 regimes in Fig.23, a e [1.69614, 1.6975], a e [1.6975, 1.698] and 
a e [1.698, 1.69847]. Typical example for each regime are shown in Fig. 3. 

To see the mechanism of the difference of dynamics in these parameter region, 
it may be useful to adopt the invariant measure Peq,ni^) the effective map, 

) = F„+6 o F„+5 o • • • o F„+i o F„, (19) 

as we have already introduced in Section 6. In Fig. 24, three examples of p^^nl^) 
are plotted as a function of time, corresponding to the three regimes mentioned 
above. In Fig. 24(a), seven distinct regions arc stabilized successively. For this 
parameter, the effective nonlinearity parameter A is close to, but smaller than, 
the tangent bifurcation point of the period seven window in the logistic map. 
Therefore if the fluctuation of the mean field were ignored, none of the seven 
regions would be stabilized because the seventh iterate of the logistic map 
Eq.(19) does not cross with y = x. With the mean field dynamics, on the 
other hand, the effective map Eq.(19) is modified to cross with y = x at 
a few regions where p^J^X^) > O(for instance between n = 2000 and 2100 
in Fig. 24 (a)). In Fig. 24 (a), two or three regions are stabilized. After some 
duration, these regions come to be destabilized again by crisis (for instance 
at n ~ 2100 in Fig. 24(a)). After the crisis Pegln(^) spreads over the whole 
region because none of the seven regions of the map Eq.(19) cross with y = x. 
Then stabilized regions switch to different positions. This process continues 
successively. 

When the parameter a is increased, the number of regions stabilized by the 
tangent bifurcation of the map Eq.(19) is increased (Fig.24(b)). In Fig. 24(b), 
5, 6, or 7 regions are stabilized successively. This corresponds to the second 
regime in Fig.23. 
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Fig. 23. (a)Mean square deviation and (b)rotation number of the mean field dy- 
namics arc plotted as a function of a. For the parameter a = 1.6962, 1.69755, and 
1.69844 indicated in the figure (a), see Fig. 3 and Fig. 24 
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Fig. 24. Invariant measure of the effective map Eq.(19) peq,n{x) is plotted as a 
function of time. The horizontal axis is time and the vertical axis is x. In this figure, 
the region with peq\n{x) > is plotted by a solid line. When the whole the region 
is filled by a line, none of the seven regions of the map Eq.(19) cross with y = x, 
and all the regions are connected as a single attracting set. In (a) and (b) some 
of the seven regions of the map Eq.(19) cross with y = x, while the other regions 
do not. In (c), on the other hand, all seven regions of the map Eq.(19) cross with 
y = X, while some of the seven regions are destabilized by crisis. The parameters 
correspond to that of Fig.3. 
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With the further increase of a, all the seven regions of the map Eq.(19) al- 
ways cross with y = x, while some of these seven regions are destabilized 
by crisis(Fig.24(c)), as we have shown in section 6. With the increase of a, 
time duration of crisis in each seven region is increased, and all the seven 
regions start to be destabilized by crisis at the same time step. Then this 
tongue structure ends up and collective dynamics in the p7 tongue structure 
is collapsed (at a = 1.69847). Then the amplitude of mean-filed dynamics is 
reduced smaller to about ]|^(see Fig. 23). 

Although we have explained the bifurcation in the internal tongue structure for 
the period-7 case, this kind of bifurcation structure is common to a band region 
in any period. For instance, in Fig. 13(a) with the period-3 tongue structure 
(starting from A e [1.75, 1.79032] at e = 0.0) and in the period-5 tongue 
structure(starting from A e [1.6244,1.6333] at e = 0.0), similar bifurcation 
structure can be observed, where the change in the number of coexisting stable 
regions makes such bifurcation structure. 



7.3 Scaling of Tongue Structures 

As we have shown in §5.2, the width in the parameter A of each tongue 
increases with e (see Fig. 13 (a)). Here, we discuss the scaling of each tongue 
structure. 

In general, the effective map for each element is given by 

Fjf^ = Fn+p-l O Fn+p-2 o • • • O o F„, (20) 



with Fn{x) — (1 — e)(l — ax'^) -\- ehn- By adopting the effective nonlinearity 
parameter A and rescaling of x„(i), which we have introduced in §5.2, F„(x) 
may be transformed into 



where 5hn = hn — {h). As we have mentioned adobe, an internal bifurcation 
condition, e.g. crisis bifurcation, and tangent bifurcation, is satisfied for some 
steps within some parameter interval 74s^ai;(e) < A{e) < Aiargei^)- (Indeed this 
region corresponds to each tongue structure). Roughly speaking, Sh^ varies in 
time within [— vMSD : vMSD] for a given parameter. Thus, the minimum 
and maximum parameter of A in a tongue structure is a function of —ey/MSD 
and -l-eVMSD respectively, i.e., Agmaii and Aiarge are given as 
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Fig. 25. Scaling relation of tongue structure for period-2 (a), period-3(b), and pe- 
riod-5(c). edh are plotted as functions of bA, where 6h = VMSD and bA indicate 
the deviation from the band merging point(^ = 1.543689012692076)(Fig.(a)), the 
crisis bifurcation point of period-3 window(A = 1.790327491999345)(Fig.(b)), and 
the crisis bifurcation point of period-5 window(A = 1.633358703691276)(Fig.(c)) of 
the logistic map, respectively. Line in each figure is proportional to bA. Hence, the 
edge of A in a tongue structure varies linearly with tbh. The width of a tongue 
structure increases proportional to tbh. 



A,raaii{e) = + A^{-e^/y^) + O ((e ■ ^MSD)^) , 
Aiarae{^) = ^0 + ^i(+e\/MSD) + O ((e • \/MSD)2) . (22) 

where Aq is a bifurcation parameter of the logistic map, and is positive 
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Fig. 26. The maximum and next maximum value of Sh = VMSD in each tongue 
structure are plotted. They are obtained by sampling the data of MSD by changing 
A in each tongue with for each e. 

constant. 



For instance, the band merging point of logistic map separates into two lines 
as is shown in Fig.22. In Fig. 25, tongue structures around the parameter A 
at the crisis bifurcation point of the period-2 window(band merging point), 
period-3 window, and period-5 window are shown. The linear scaling relation 
with evMSD is clearly seen as to the change of Asmaii and Aiarge- Thus the 
width of the tongue structure grows linearly with e • Sh. 

To obtain the scaling of the tongue structure as a function of e, we have to 
know the dependence of MSD on e. In Fig. 26, the growth of square root of 
the MSD of the mean field dynamics in a tongue structure with the coupling 
strength e is shown for several tongue structures. While it has been pointed 
out that the amplitude of the mean field dynamics may grow linearly with 
e for globally coupled logistic map [8, 19], Fig. 26 indicates a possibility that 
there is a deviation from the linear scaling with e for the amplitude of the 
mean field. It should be noted that while we have payed attention mainly to 
tongue structures relevant to windows of the logistic map, in the reference[19] 
such window structures in the logistic map are out of consideration. In other 
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Fig. 27. Time scries of elements. Elements are accumulated to three bands. In the 
figure time series of 100 elements out of 10^ are plotted at every three steps. A lot 
of attractors are realized depending on the population ratio to each band, while 
elements are desynchronized each other. The parameters are a = 1.85, e = 0.018, 

= 10^ 

words, they have focused on the collective dynamics arising from completely 
chaotic dynamics in the logistic map. Note also that our analysis is based on 
the rescaled nonlinearity parameter A, while the studies in the references [8,19] 
are based on a. Possible distinction between the collective motions originated 
in chaos and window will be discussed in §9 again. 



8 Hysteresis, Multiple Attractors, and Coexistence of Different 
Types of Motion 



Even if the control parameters are same, depending on its initial condition, 
there can exist more than one attractors of the collective motion. 

Most straightforward examples of multiple attractors are given with the use 
of band splitting. At same parameter region, while there is no mutual syn- 
chronization, elements are accumulated to few bands and never change their 
band(Fig.27). Therefore, a lot of attractors are realized depending on the 
population ratio to each band, as long as stability conditions are satisfied. For 
instance. In Fig. 27, elements are accumulated into three bands, (see e.g. [14] 
for the case with a tent map) . 

The next example of multiple attractors is concerned with hysteresis phenom- 
ena of collective motion, which can be observed at the edge of the tongue 
structure in the parameter space. In Fig. 28, hysteresis curve of MSD is ob- 
served by increasing or decreasing the control parameter a with the use of 
the final state of a simulation at the previous value of a as the next initial 
condition. Thus in a e [1.69848 : 1.69858], at least two different attractors of 
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Fig. 28. Hysteresis curve observed by increasing or decreasing progressively the 
control parameter a while keeping the final state of a simulation at given a as the 
initial condition for the neighboring value a — 6a{x) and a + 5a{+). e = 0.008. The 
MSD calculated starting from a random initial condition are also plotted (o). 

the collective motion coexist. In Fig. 29, time series and return map for each 
attr actor are shown. Note that, in this case there is no separated bands in 
contrast with the previous examples. 

The third example of multiple attractors, one attractor has a band structure 
(Fig. 27) and the other not (Fig.30). For the former attractor, elements are 
accumulated in a few bands, while for the latter elements spread over whole 
range of x. Moreover for the former type, there exists a lot of attractors with 
a different ratio of population in each band, as in the first example. 

Another important topic related to the multiple attractor is the coexistence 
of different kinds of element motions. When the elements arc accumulated 
into few bands, depending on the ratio of population in bands, the motion of 
elements in each band is different. In Fig.31, two kinds of element motions arc 
plotted for attractors with different population splitting ratio into bands. At 
this parameter, there is a three-band structure, and elements are accumulated 
into two groups of these three bands. Note that while these groups are similar 
to clusters (cf. [6]), but the value of elements Xn{i) in each group are different 
each other. Depending on the ratio of population in the two groups, two kinds 
of element motion coexist. Relevance of such coexistence to the problem of 
cell differentiation is discussed in [31]. 
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Fig. 29. Time series and return map for two attractors in Fig.28. The parameters 
are a = 1.69855, e = 0.008. 

9 Summeiry and Discussions 



In the present paper, we have studied the collective motion in desynchronized 
state of globally coupled logistic maps. It is shown that the motion with a much 
longer time scale and lower dimension can emerge in macroscopic dynamics, 
such as the mean field dynamics. The amplitude of collective motion (mean 
square deviation of the mean field distribution) is studied by changing the 
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Fig. 30. Time series of elements for the same parameter in Fig.27 starting from 
different initial condition. In contrast with Fig.27, elements spread over x. The 
mean field dynamics for this time series shows quasi-periodic-like motion. Time 
series are plotted for 100 elements at every three steps. The parameters are a = 1.85, 
e = 0.018, N = 10^. 
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Fig. 31. Return maps of element motions to show the coexistence of different kinds 
of motions. Pi (circled) and P2 represent the element motion belonging to a dif- 
ferent band in the dynamics, respectively. Two kinds of motions (Pi, and P2) co- 
exist. Pi : P2 in each figure indicates the number of population in each group. 
With the change of ratio, two kinds of motions are varied. The parameters are 
a = 1.88, e = 0.04, = 10^. 
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Fig. 32. Schematic diagram of tongue bifurcation structure. Regions that grow from 
a bifurcation point of the logistic map constitute a tongue structure. See §7 

nonlinearity parameter a and the coupling strength e. By introducing the 
effective nonlinearity parameter A with rescaling of Xn{i), tongue structure 
has been detected in {A, e)-plain. Each tongue structure corresponds to a 
periodic window of logistic map. 

Focusing on the tongue structure, we have demonstrated how such a collec- 
tive motion emerges. Self-consistent dynamics between the mean field dynam- 
ics(macroscopic dynamics) and each element (microscopic dynamics) is found 
to be formed, so that such a collective motion is possible. This self-consistent 
dynamics is formed by the following circulation: accumulation of elements 
into some regions leading to the change in the mean field dynamics, which 
introduce the stabihty change of the the regions, and flow of elements into 

a different region, which, again This gives internal bifurcation in elements 

and in time. 

The bifurcation is also seen in the parameter space. Since the nature of the 
internal bifurcation varies with the nonlinearity parameter a in a tongue struc- 
ture, the number of coexisting regions in x changes, which makes the collective 
motion qualitatively different. Hence, in a tongue structure, different kinds of 
collective motions have been observed. (A schematic figure of tongue structure 
is presented in Fig. 32). 

With the increase of the coupling strength e, each tongue structure grows in 
proportion to e ■ Sh, where dh is the amplitude of the mean field variation. 
Hence the width of each tongue increases with if 6h oc e. In contrast with 
earlier studies[8,19] supporting this linear scaling, however, our calculation 
suggests that the scaling may obey a different power law. 
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Existence of multiple attractors with a different collective motion have also 
been reported. This means that there can be non unique self-consistent dy- 
namics between microscopic and macroscopic motion. It should be noted that 
in the case of the multiple attractors, since the average mean field is clearly 
different by attractors, they take different values of the effective nonlinearity 
parameter A, and are distinguishable clearly in the {A, e)-plain. 

Our tongue structure is based on the windows that exist in the logistic map. 
Since the window exists in any neighborhood in the parameter space of the 
logistic map and the width of the tongue structure increases with the coupling 
e, the tongue structure is expected to occupy a relatively large region in the 
parameter space. This is one of the reasons why we have focused our atten- 
tion on the collective behavior in the tongue structures. Still, we have to note 
that there is a positive measure in the parameter space of the logistic map, 
corresponding to chaos. Hence, at least at small coupling regime in our GCM, 
there are parameters with a positive measure which do not belong to any 
tongue structure. Indeed, we have observed that the amplitude of the mean- 
filed variation drops less than to l/10(see Fig. 12), at the parameter where the 
tongues structure disappears. Although no clear structure in the return map 
is detected there, this motion again has some hidden coherence and is distin- 
guishable from noise. Analytical estimate of the mean field dynamics by S. V. 
Ershov, et al.[19] is expected to correspond to such chaos-originated regime. 
However, we need further study to clarify the mechanism of the collective 
dynamics there, and characterize the high-dimensional chaos. 

While in the desynchronized state elements are completely desynchronized 
each other and all the Lyapunov exponents are positive, for some parameter 
regime, some kind of predictability may emerge in the macroscopic variables. 
However, since low-dimensional (0(1)) collective dynamics has not been ob- 
served, microscopic and macroscopic dynamics are not separated completely. 
This is why we need a self-consistent description between microscopic variable 
(x„(z)) and the mean- field. On the other hand, it might be also important 
to study how such characteristics of the collective motion are reflected on A^- 
dimensional phase space structure, or on microscopic quantities, such as the 
Lyapunov spectrum. With such study, the mechanism for our collective mo- 
tion must be clearly distinguishable from the self-organization mechanism [32] 
or the slaving principle [33]. Although we have presented a heuristic way to 
extract such self-consistent dynamics in the present paper, it is hoped that 
a systematic method to characterize the (high-dimensional) collective motion 
will be developed in future^. 



So far, we have no conventional tool for detecting the lower dimensional collective 
signals out of high dimensional signals. In [34], we will develop a tool to distinguish 
and characterize several collective dynamics in GCM. 
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